Exact critical points of the 0{n) loop model on the martini and the 3-12 lattices 
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We derive the exact critical line of the 0(n) loop model on the martini lattice as a function of the 
loop weight n. A finite-size scaling analysis based on transfer matrix calculations is also performed. 
The numerical results coincide with the theoretical predictions with an accuracy up to 9 decimal 
places. In the limit n — > 0, this gives the exact connective constant fi = 1.7505645579 ■ • ■ of self- 
avoiding walks on the martini lattice. Using similar numerical methods, we also study the 0(n) 
loop model on the 3-12 lattice. We obtain similarly precise agreement with the exact critical points 
given by Batchelor [J. Stat. Phys. 92, 1203 (1998)]. 
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Introduction. The 0(n) loop model [H originates from 
the high-temperature expansion of the 0(7i) spin model 
. It can be considered a model describing a noninter- 
secting loop gas. On lattices with coordination number 
three, the partition function is very simple: 



Z : 
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where the sum is over all configurations of non- 
intersecting loops denoted as Q; x is the weight of a 
bond (an edge occupied by loop segments), or an occu- 
pied vertex, and n the weight of a loop, b is the number 
of bonds or occupied vertices, and I the number of loops. 

Generally speaking, there is a high-temperature phase 
with dilute loops and a low-temperature phase with 
dense loops. At the transition point Xc, the longest loop 
grows to infinity and begins to percolate the system. 
The critical properties of this transition are universal, 
which are well described by the Coulomb gas theory 
0. However, the determination of the critical points of 
this model on various lattices remains to be treated case 
by case. Exact 0(n) critical lines have been found on 
the honeycomb lattice l^-Q ; the square lattice [7] and 
the triangular lattice [8|. In addition to this transition 
point, there are several other branches of critical be- 
havior, e.g., 'branch 0', which describes a higher critical 
point, as reported in Refs. [8l-[l0|. 

In the n — >■ limit, the critical 0(n) loop model de- 
scribes long polymers in a good solvent or self-avoiding 
walks (SAWs)[llj. The study of the 0(n) loop model 
has led to a wealth of information on the configuration 
properties of S AWs [l2| • The number of configurations 
of SAWs in k steps, i.e., Ck scales as [ll| 



(2) 



for large k. j — 43/32 is a universal critical exponent, 
which can be obtained via the Coulomb gas theory Q- 
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/i is the connective constant which is lattice dependent, 
and equals to 1/xc of the n — > loop model. Although 
the studies on the SAWs have advanced a lot since it 
was introduced [l3l, the values of /i for most of two- 
dimensional lattices are found numerically ,15,- 19] . But 
the exactly known critical line of the 0(n) loop model 
on the honeycomb lattice 3 provides the exact result 

/i = v2 + \/2 for the model on the honeycomb lattice. 
The critical line of the honeycomb 0(n) loop model was 
found by an exact mapping on a Potts model [ij and by 
using the Bethe Ansatz 5, 6]: 
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The phase diagram inferred from this result has been 
well verified by different numerical methods [20l. [2ll |. 
Batchelor derived the critical points of the 0(n) loop 
model on the 3-12 lattice by mapping the honeycomb 
loop model to the 3-12 lattice j23|: 



1 -I- x^n 
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This result is in agreement with existing exact results 
[22I I for the cases n = 1 (Ising) and n — Q. Making 
use of this result, Batchelor obtained /i — 1.711041 • • • 
for the 3-12 lattice, which coincides with the result pre- 
viously found by Jensen and Guttmannfisj using other 
methods. In the present paper we shall provide some in- 
dependent numerical results for the critical point of the 
0(n) loop model on the 3-12 lattice, and on its phase 
behavior. 

Inspired by Batchelor's work, we studied the 0(n) 
loop model on the martini lattice. This lattice was 
first proposed by ScuUardjll] in the study of percola- 
tion. The percolation threshold and the critical points 
of the g-state Potts modeljl^, HH on this lattice are 
known exactly |2^], but the critical points of the 0(n) 
loop model are not known yet. In this paper, we derive 
the exact critical points of the 0(n) loop model as a 
function of n on the martini lattice. In the limit n — > 0, 
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the result gives the exact connective constant of SAWs 
on the martini lattice. In addition, we build the transfer 
matrix (TM) and apply a finite-size scaling analysis for 
a numerical study of the model on the martini and the 
3-12 lattice. 

Critical points of 0(n) loop model on the martini lat- 
tice. Consider a honeycomb lattice with two sublattices 




(a) 



(b) 



FIG. 1. (Color online) The 0(n) loop model on the honey- 
comb lattice (a) and the martini lattice (b). 




FIG. 2. (Color online) Mapping of the vertex configurations. 
The weights of the different vertices are shown. The same 
weights apply to versions rotated by ±27r/3. 

A and B. A loop configuration Q also denotes the oc- 
cupations of sublattice A, as shown in Fig. [IJa). We 
rewrite the partition function ([T|) in the following way: 



(5) 



where v = b is the number of vertices visited by loop 
segments, Va and va — v/2 are the number of vertices 
and the number of visited vertices of sublattice A, re- 
spectively. Thus the weight of a visited A vertex shown 
in Fig. IHa) is x^, the weight of an empty A vertex 
shown in Fig. [Hd) is 1. 

Now consider the 0(n) loop model on the martini lat- 
tice, which is constructed by replacing the 'star' (an A 
vertex) shown in Fig. O^d) by the structure shown in 
Fig. m^e). Each occupied (empty) A vertex corresponds 
to two possible occupations on that structure, as shown 
in Fig. [2] (b) and (c) ((e) and (f)). Thus, any given 
configuration Q of loops on the honeycomb lattice maps 
to the sum of 2^-* possible loop configurations on the 
martini lattice. Let x be the weight of a bond for the 
martini loop model, the partition function of the loop 
model on the martini lattice can be obtained by sum- 
ming on Q: 



^ 1 x^n 

G 



(6) 



Mapping {x'^ + x^)/{l -\- x^n) x^, we obtain in 
([S]) multiplied by a trivial factor. It follows the critical 



points Xc of the 0{n) loop model on the martini lattice 



(7) 
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This result is in agreement with an existing result for 
the 0(1) loop model, which is equivalent to the high 
temperature expansion of the Ising model model j27| 
with X = tanhX^, where is the coupling of Ising 
spins sitting on the vertices of the martini lattice. Ac- 
cording to (0), K[ = 0.749790959036 ••• , which coin- 
cides with the exactly known critical point of the q = 2 
Potts model on the martini lattice [26| . 

As another special case, we present the exact connec- 
tive constant of the SAW on the martini lattice. The 
substitution /i = l/xc in ([7]) for n = determines /x as 
the solution of 



1 1 _ V2 
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which yields n = 1.750564567897 • • • . 

We may further generalize above results by allowing 
bonds on the small triangles in Fig. ^e) have weight 
(xt) different from those on the remaining ones (xg). 
Following the mapping described above, we thus obtain 
a critical line in the xt versus Xg plane for a given n: 



1 -\- Xt^n 
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For the 3-12 lattice, the critical line of the generalized 
model is 



Xsjxt + xD 

1 + x^n 
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Finite-size scaling and transfer matrix calculation. 
Consider the lattice (the 3-12 or the martini lattice) 
wrapped on a cylinder with circumference L. The mag- 
netic correlation function of the 0(n) spin model is 
translated as the probability that two sites at a distance 
r are linked by a single loop segment 0], gr = Z' /Z, 
where Z' = Y^g, x^n\ and Q' denotes the configura- 
tions that connect sites and r by precisely one single 
loop segment. In our transfer-matrix analysis of the 
finite-size-scaling behavior, it is sufficient to substitute 
the configurations Q' that connect any site of row to 
any site of a row at a distance r as measured in the 
length direction of the cylinder. The exponential decay 
of gr at large distances is determined by the magnetic 
gap in the eigenvalue spectrum of the transfer matrix. 
The scaled magnetic gap Xh{x,L) = || ln(A(")/A(i)), 
where A^°) and A'-'^' is the largest eigenvalue of the TM 
for Z and Z' , respectively. C is a geometrical factor de- 
termined by the ratio between the unit of L and the 
thickness of a row added by the transfer matrix. An- 
other scaled gap Xt{x,L) = ^ln{A^°'> / A^'^^), describes 
the exponential decay of the energy-energy correlation, 
with A^^^ the second eigenvalue of the TM for Z. 

The TM techniques of the 0(n) loop model are well 
described in the literature, e.g., see [§]. The procedure 
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of sparse matrix decomposition for the martini and the 
3-12 lattice equals to that for the honeycomb lattice 
with the addin g un its suitably chosen [20|. For further 
details, see [li,!!]. 

According to the finite-size scaling ^3d\ and the con- 
formal invariance 31] theory, the scaled gap Xi{x,L), 
in the vicinity of the critical point, satisfies 

X^{x, L)^X^+ a{x - Xc)i^' + buL^^ + • • • , (11) 

where Xi{i ~ h,t) is the magnetic and the temper- 
ature scaling dimension, respectively; yt is the ther- 
mal exponent; u denotes the leading irrelevant field, 
Uu is the associated irrelevant exponent, a, b are un- 
known constants. Such behavior is illustrated in Fig. 
[3]^ a) and (c) for the case n = 0. The critical point 
is estimated by numerically solving x in the equation 
Xh{x, L) = Xh{x, L—1), with system sizes up to L = 16. 
The solution Xc{L) converges to the critical value Xc as 

Xc{L) = + a'uL^'^-y' + ■■■ , (12) 

where a' is an unknown constant. The numerical esti- 
mations and the theoretical predictions of the critical 
points of the martini lattice are listed in Table HI Our 
numerical estimations coincide with the theoretical pre- 
dictions in 9 decimal numbers. 

The universal values of Xi and the conformal anomaly 
c of the two-dimensional 0(n) loop model are exactly 
known as [1, [3l| 



X.^i^, ,13) 

2.9 8 g 

where n = — 2cos(7r5), 1 < 5 < 2. 

At criticality, Xi{L) converges as follows to Xi with 
increasing system size L 

X,{L)^X, + b'Ly- + -- - , (14) 

where b' is an unknown constant. The free energy den- 
sity /(L) = ClnA^^V^ scales as [Hill] 

f{L) = /(oo) + ^. (15) 

We then calculate f{L), Xt{L) and Xh{L) for a se- 
quence of systems up to size i = 16 at the critical points 
d?]). Fitting the data according to (jH)) and ((T5|) . we ob- 
tain the scaling dimensions Xh^ Xt and the conformal 
anomaly c, which are also listed in Table HI Our numer- 
ical estimations are in agreement with the theoretical 
predictions with a high accuracy. For Xh , Xt , our re- 
sults are also consistent with the Monte Carlo results 
forn > 1 ill,!!!. 

When n approaches 2, yt — 2 — Xt 0. The cor- 
rections to scaling due to the leading irrelevant field 
become relatively strong. At n = 2, yt is exactly 0, 
so that intersecting points between the curves Xi{x,L) 



and Xi{x,L — 1) may be absent, as already suggested 
by pT|) . and as indeed observed in Fig. [3][b) and (d). 
Therefore, we can't numerically determine the critical 
point Xc in the usual way. However, c and Xh,Xt are 
still estimated at the theoretical critical point. 

Similar analysis is also performed to the 3-12 lat- 
tice. Theoretical predictions and numerical estimations 
of critical points for several values of n, which agree in 
a high accuracy, are listed in Table [TTl ^From Tables |T] 
andim we can see that the values of Xh, Xt and c for a 
two-dimensional 0(n) loop model on the martini lattice 
coincide with those of the 3-12 lattice, as expected by 
the hypothesis of universality. 
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FIG. 3. (Color online) Scaled gaps Xh{x, L), Xt{x, L) versus 
X for the 0(n) loop model on the martini lattice with L=8 
to 16. (a) and (c): n — 0; (b) and (d): n = 2, with Xc = 
0.840896 indicated. Lines connecting data points are added 
to guide the eye. 

Conclusion. We derived the exact critical points of 
the 0(n) loop model on the martini lattice, and per- 
formed a finite-size scaling analysis based on numerical 
TM calculations. Our numerical estimations agree with 
the theoretical predictions, within a margin that is typ- 
ically of order 10~^. This rather high precision may be 
related to the vanishing of the leading irrelevant field in 
the Nienhuis result Q for the critical line. 

In the limit n — > 0, the critical point gives the exact 
connective constant /i = 1.7505.... of the SAWs on the 
martini lattice. 

The exact critical points of the 0(n) loop model on 
the 3-12 lattice derived by Batchelor are also verified. 

The conformal anomaly, the magnetic and the tem- 
perature scaling dimensions of the 0(n) models on the 
two lattices are numerically calculated. The estimations 
coincide with the theoretical predictions, as expected 
according to the universality hypothesis. 
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TABLE I. Critical points Xc, conformal anomaly c, magnetic and temperature scaling dimensions Xh,Xt of the two- 
dimensionaJ 0(n) loop model on the martini lattice. (T=Theoretical prediction, N=Numerical estimation.) 



n 


j-cv J- ; 












V J- 7 







0.571244285 


0.571244285(1) 








0.1041667 


0.104166(1) 


2/3 


0.666668(2) 


0.25 


0.584248605 


0.584248605(1) 


0.1300704 


0.1300705(2) 


0.1100192 


0.1100193(1) 


0.7395254 


0.739526(1) 


0.5 


0.598867666 


0.59886768(2) 


0.2559499 


0.255950(1) 


0.1154420 


0.1154420(1) 


0.8177559 


0.817756(1) 


0.75 


0.615559079 


0.6155590(1) 


0.3788781 


0.3788783(3) 


0.1204452 


0.1204452(1) 


0.9035105 


0.9035104(2) 


1.0 


0.635024224 


0.63502422(1) 


0.5 


0.500000(1) 


0.125 


0.12500000(1) 


1 


1.00000(1) 


1.25 


0.658437850 


0.65843786(1) 


0.6205051 


0.6205053(2) 


0.1290128 


0.1290127(1) 


1.1126008 


1.1126007(1) 


1.5 


0.688067393 


0.68806739(2) 


0.7418425 


0.741842(1) 


0.1322435 


0.1322434(1) 


1.2518912 


1.25189(1) 


1.75 


0.729662053 


0.729664(3) 


0.8662562 


0.8662563(1) 


0.1339623 


0.13396(1) 


1.4457176 


1.445718(1) 


2.0 


0.840896415 




1 


1.000000(1) 


0.125 


0.1250000(1) 


2 


2.00000(1) 



TABLE IL Critical points Xc, conformal anomaly c, magnetic and temperature scaling dimensions Xh,Xt of the two- 
dimensional 0(n) loop model on the 3-12 lattice. (T=Theoretical prediction, N=Numerical estimation.) 



n 




a^c(N) 


c(T) 


c(N) 




Xh{N) 




Xt(N) 





0.584439429 


0.584439429(1) 








0.1041667 


0.104167(1) 


2/3 


0.666668(1) 


0.25 


0.601034092 


0.601034092(1) 


0.1300704 


0.1300705(2) 


0.1100192 


0.1100193(1) 


0.7395254 


0.739526(1) 


0.5 


0.620240607 


0.62024060(1) 


0.2559499 


0.255950(3) 


0.1154420 


0.115442(1) 


0.8177559 


0.817756(1) 


0.75 


0.642967899 


0.6429678(1) 


0.3788781 


0.3788783(2) 


0.1204452 


0.1204452(1) 


0.9035105 


0.903510 


1.0 


0.670697664 


0.67069766(1) 


0.5 


0.499999(1) 


0.125 


0.1250000(1) 


1 


1.000000(1) 


1.25 


0.706102901 


0.70610291(1) 


0.6205051 


0.620505(1) 


0.1290128 


0.1290127(1) 


1.1126008 


1.112600(1) 


1.5 


0.754845016 


0.754845(1) 


0.7418425 


0.741842(1) 


0.1322435 


0.1322435(1) 


1.2518912 


1.25189(1) 


1.75 


0.833205232 


0.83320(2) 


0.8662562 


0.86626(1) 


0.1339623 


0.13396(1) 


1.4457176 


1.445718(2) 


2.0 


l.i72r):5i()77 




i 


0.91)1)1)99(1) 


0.12") 


0.12.")()()()()(1) 




2.00000(1) 
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